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We propose the use of optical lattice clocks operated with fermionic alkaline-earth-atoms to study 
spin-orbit coupling (SOC) in interacting many-body systems. The SOC emerges naturally during 
the clock interrogation when atoms are allowed to tunnel and accumulate a phase set by the ratio of 
the “magic” lattice wavelength to the clock transition wavelength. We demonstrate how standard 
protocols such as Rabi and Ramsey spectroscopy, that take advantage of the sub-Hertz resolution 
of state-of-the-art clock lasers, can perform momentum-resolved band tomography and determine 
SOC-induced s-wave collisions in nuclear spin polarized fermions. By adding a second counter- 
propagating clock beam a sliding superlattice can be implemented and used for controlled atom 
transport and as a probe of p and s-wave interactions. The proposed spectroscopic probes provide 
clean and well-resolved signatures at current clock operating temperatures. 


The recent implementation of synthetic gauge fields 
and spin-orbit coupling (SOC) in neutral atomic gases 
m is a groundbreaking step towards the ultimate goal 
of using these fully-controllable systems to synthesize and 
probe novel topological states of matter. So far optical 
Raman transitions have been used to couple different in¬ 
ternal (e.g. hyperfine) states while transferring net mo¬ 
mentum to the atoms. However, in alkali atoms Raman- 
induced spin-flips inevitably suffer from heating mech¬ 
anisms associated with spontaneous emission. While 
this issue has not yet been an impediment for the in¬ 
vestigation of non-interacting processes or mean field ef¬ 
fects H E], it could limit the ability to observe inter¬ 
acting many-body phenomena that manifest at longer 
timescales. Finding alternative, more resilient methods 
for generating synthetic SOC, and probing its interplay 
with interactions, is thus highly desirable. 

To reduce heating, the use of atoms with richer inter¬ 
nal structure such as alkaline-earth-atoms (AEAs) [3l |7] 
or lanthanide atoms such as Dy and Er [8] have been 
suggested. Here, we propose an implementation of SOC 
using cold AEAs in an optical lattice clock (OLC) [9] and 
discuss how SOC physics can be probed in interacting 
many-body systems. Our proposal relies on the fact that 
SOC naturally arises in OLCs because the clock laser 
imprints a phase that varies significantly from one lat¬ 
tice site to the next as it drives an ultra-narrow optical 
transition. Our proposal offers several advantages over 
prior schemes. It uses a direct transition to a long-lived 
electronic clock state with natural lifetime > 10^ s [TO], 
and thus heating from spontaneous emission is negligi¬ 
ble. Second, our proposal takes advantage of the sub-Hz 
resolution of clock lasers Einally, we can probe 

the interplay of interactions and SOC by operating in 
the regime where the interaction energy per particle, U, 
is weak compared to the characteristic trapping energies 
[T 5 H 2 O] but comparable to SOC scales determined by J, 
the tunneling, and the clock Rabi frequency. 

We propose three protocols for probing SOC in an 
OLC and demonstrate that all provide clear signatures 
under current operating conditions. In the first proto- 



FIG. 1: (Color online) (a) A clock laser along the Z direction 
of wavelength A and Rabi frequency Q interrogates the ^So{g)- 
^Po(e) transition in fermionic alkaline-earth atoms trapped in 
an optical lattice with magic wavelength Am- The transverse 
confinement is provided by Gaussian curvature of the lattice 
beams with harmonic frequency hiy ±. Many transverse modes 
n are populated at current operating temperatures, (b) The 
phase difference <f between adjacent sites j and (j + l) induces 
SOC when atoms can tunnel with mode-dependent tunnel¬ 
coupling Jn, realizing a synthetic two-leg ladder with flux (f> 
per plaquette. 


col, we demonstrate that momentum-resolved tomogra¬ 
phy of chiral band structures can be performed using 
Rabi spectroscopy in the parameter regime J ^ ^ I/. 

In the second protocol we show that the modification 
of collisional properties by SOC [21] is manifest in 
standard Ramsey spectroscopy, focusing on the param¬ 
eter regime J > I/. In the final protocol, controlled 
and spatially resolved atomic transport [231130] is in¬ 
duced by an additional counter-propagating clock beam 
that exhibits a controllable phase difference with respect 
to the original probe beam, and interaction effects be¬ 
yond mean field modify the dynamics in the parameter 
regime Cl ^ J ^ U. Although inelastic collisions in 
the excited state impose limitations in the probing time 
[Ii|20l|3ll[34], we also show that they can be used as a 
resource for state preparation and readout. 

SOC implementation Current OLCs interrogate the 
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^So{g)-^Po{e) transition of ensembles of thousands of nu¬ 
clear spin polarized fermionic AEAs trapped in a deep 
ID optical lattice that splits the gas in arrays of 2D pan¬ 
cakes [9] (see Fig 0 . The lattice potential uses the magic- 
wavelength, to generate identical trapping conditions 
for the two states. At current operating temperatures, 
T ^ juK [15] the population of higher axial bands is negli¬ 
gible (< 5%). On the other hand, along the transverse di¬ 
rections, where the confinement is provided wholly by the 
Gaussian curvature of the optical lattice beams, modes 
are thermally populated with an average number of mode 
quanta (n) ^50. To generate SOC coherent tunneling 
between lattice sites is required. Our proposal is to su¬ 
perimpose a running-wave beam on the lattice potential: 

Vonst + ^COS2(^)] ; 

( 1 ) 



FIG. 2: (Color online) (a) SOC band structure for 6 = —2Jo, 
Do = Jo (solid lines) and Dq — 0-05Jo (dashed lines). The ax¬ 
ial depth is VIEr = 12, Jo = 42 Hz and ~ 900 Hz. Col¬ 
ors correspond to state character, with g (e) being more blue 
(red), (b) Chiral Bloch vector angle, Ooq, in the xz plane ex¬ 
tracted from Rabi spectroscopy using the protocol explained 
in the text [38]. The figure shows three temperatures for the 
parameters of (a). 



this increases the transverse confinement without signif¬ 
icantly affecting the axial motion [35]. Here, wq is the 
beam waist, R the transverse radial coordinate, Z the ax¬ 
ial coordinate, V the axial lattice corrugation, and Konst 
the running-wave induced potential. By increasing Konst 
as V is lowered, the transverse confinement frequency 

iy± ~ '\J (Vonst + V) is kept constant while the 

tunneling rate along the axial direction increases. Here 
m is the atom mass. 

Since Kxt (R) is discretely translationally invariant 
along the axial direction, atoms trapped in the low¬ 
est axial lattice band are governed by the Hamiltonian 

= m 


q,a 




E 


Dn 


-hH.c. 


.( 2 ) 


Here and throughout q = qa is the dimensionless product 
of axial quasimomentum q and lattice spacing a = 

^a,n,q annihilates a fermion in the two-dimensional trans¬ 
verse mode n, quasimomentum q, state a = ±(for e and 
g), and The energy {q) = 

q;| + — 2 J„cos(g) has contributions from the mode 

dependent tunneling Jn, the average energy of the trans¬ 
verse mode n, E^, and the laser detuning d. Dn is the 
Rabi frequency for mode n. The clock laser with wave¬ 
length A imprints a phase that varies between adjacent 
lattice sites by 0 = ttA^/A. 

If ones views the two internal states as a discrete syn¬ 
thetic dimension [saiszj, as shown in Fig.[^ describes 
the motion of a charged particle on a two-leg ladder in a 
magnetic field with flux 0 per lattice plaquette. hence 
has the interpretation of many copies of those ladders - 
one for each transverse mode. By performing a gauge 
transformation ^+,g,n 5 becomes diago¬ 

nal in momentum space with the excited state dispersion 
shifted \yy (j), q ^ q (j). The latter can be then conve¬ 
niently written in terms of spin- 1/2 operators acting on 


the populated modes 

= (3) 

n,g 


where B„q = 0 , AEn ( 9 , (/>) + <5), AEn{q,(j)) = 

2 J„ [cos (g) — cos (g + (/))], and are spin- 1/2 an- 

gular momentum operators. The eigenstates are de¬ 
scribed by Bloch vectors pointing in the xz plane, 
with a direction specified by a single angle O^q = 

arctan ^ dependence of this angle is a 

manifestation of chiral spin-momentum locking, which is 
directly connected to the topological chiral edge modes 
of the two-dimensional Harper-Hofstadter model mm- 
An example of the eigenstates and their dispersion for 
fixed transverse mode index is given in Fig. |^a). 

Rabi spectroscopy: We now describe a spectroscopic 
protocol to probe the non-interacting chiral band struc¬ 
ture of Eq. In the limit of very weak Rabi frequency. 
Do To, the band structure is almost the “bare” one 
with D = 0 , up to a displacement of the excited state 
band by (j). Because of the displacement, there are spe¬ 
cial quasimomentum points g* where the two dispersions 
cross, see Eig. [^a). These exist for a finite window of 6 
and are signaled in the carrier linewidth. The width of 


the window is 8 Jn 



and thus when resolved it can 


be used to determine (j). At finite temperature, many 
transverse modes are populated and hence the depen¬ 
dence of Dn and Jn on n could broaden the line and in 
general prevent an accurate determination of (j). How¬ 
ever, a direct simulation of the Rabi lineshape using the 
potential Eq. Q demonstrates that the features of the 
ideal, zero-temperature lineshape are captured even for 
a temperature of 3/iK [38] . 

The ability to resolve the g* resonances that appear 
across the entire Brillouin zone (BZ) as the detuning is 
varied can be used to perform momentum resolved spec¬ 
troscopy and to precisely determine the chiral Bloch vec¬ 
tor angle Ooq for given values of Dq and 5 from Rabi 
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FIG. 3: (Color online) A clock laser pulse with pulse area 
0i imprints a phase difference cj) between atoms in neighbor¬ 
ing sites. Atom tunneling, Jn, allows for s-wave interactions, 
oc U~g^ which are signaled as a density shift in Ramsey spec¬ 
troscopy after a second pulse of area O 2 is applied. 


oscillations. We propose the use of a three spectroscopic 
sequences [38]. One sequence selectively excites atoms 
at g* from g to e and induces Rabi oscillations. An¬ 
other filters the dynamics of the excited atoms from the 
remaining g atoms and the third is used to “correct” im¬ 
perfections arising from finite temperature. As shown 
in Fig. Ilfb), this protocol allows us to extract Ooq over 
the BZ. This method is not restricted to OLCs, requiring 
only a stable probe, and complements other techniques 
for measuring band structures [401442] used in degenerate 
Fermi gases. 

Ramsey spectroscopy The standard Ramsey protocol 
starts with all atoms in the g state and subjects them 
to two strong, near-resonant laser pulses separated by a 
free evolution “dark” time r. The first pulse rotates the 
Bloch vector by an angle 0i set by the pulse area, and 
the second converts the accumulated phase during the 
dark time into a g — e population difference measured as 
Ramsey fringes. Interactions induce a density-dependent 
frequency shift in the fringes. 

Interactions between two nuclear spin polarized atoms 
depend on the motional and electronic degrees of free¬ 
dom [18]. When the atoms collide they experience 5 - 
wave interactions, characterized by the elastic scatter¬ 
ing length Oeg, when their electronic state is antisym¬ 
metric {\eg) — \g^)j\f2. They can also collide via p- 
wave interactions, described by the corresponding p- 
wave elastic scattering volumes 6^^, 6^^, and 6^^, in 
the three possible symmetric electronic configurations 
1^^), |ee), (je^f) + \ge))l\/2^ respectively. In addition to 
elastic interactions, atoms can also exhibit inelastic colli¬ 
sions. In ^^Sr only the ee type has been observed to give 
rise to measurable losses m while in yb, both ee 

and eg losses have been reported [T9[ |20l |3T] . The effect 
of losses can be compensated by tracking the population 
decay during the dark time [isms]. 

When tunneling is suppressed the differential phase im¬ 
parted by the laser is irrelevant, and as long as is the 
same for all modes - a condition well-satisfied in current 
OLCs [15] - the collective spin of the atoms within each 
lattice site remains fully symmetric after the pulse and 
only p-wave collisions occur during the dark time. Mea¬ 
surements under this condition [151 ES] indeed observed 
a frequency shift linearly dependent on the excitation 
fraction of atoms, (1 — cos 0i)/2 and fully consistent with 


a p-wave interacting model m- If instead tunneling is 
allowed during the dark time, atoms become sensitive to 
the spatially inhomogeneous spin rotation from the site- 
dependent laser phase, which in turn allows for s-wave 
collisions after a tunneling event (see Fig. |^. The s- 
wave collisions in SOC-coupled spin polarized fermions 
can lead to exotic phases of matter including topological 
quantum liquids 

In typical OLCs the interaction energy per particle 
(at the Hz level) is weaker than the vibrational en¬ 
ergy spacing, and large anharmonicity prevents resonant 
interaction-induced mode transverse changes, and so a 
description in terms of fixed transverse motional levels 
is justified [18]. In the regime of weak interactions com¬ 
pared to tunneling, we compute the dynamics perturba- 
tively, and find that SOC manifests itself in the density 
shift at short times as [38] 


Au = Az/q 


4( cos 0i sin^ | 

3h^C-xcos0i) 


(4) 


where Az/q = N{C — xcos6>i) is the density shift in the 
absence of tunneling mill], N the mean atom num¬ 
ber per pancake, {J^)tr the thermally averaged squared 
tunneling rate, ( = /2, x = — 

2H^^)/2, and C = (W^-H^^)/2, with = bl^,{P)TR 

and = ae^(5')TK-Here {P)tr ^ and {S)tr cc 

correspond to the thermal averages of the p-wave 
and s-wave mode overlap coefficients respectively [18], 
and Tr is the radial temperature. For the JILA ^^Sr 
clock operated at Tr ^ (1 — 5)/iK and 6>i ^ tt for r ^ 80 
ms, Az/q ^ —5Hz. Here, Er is the recoil energy. Since 
SOC introduces contributions from s-wave interactions, 
which can be one order of magnitude larger than p—wave 
at Tr ^ 1/iK, then from Eq. [^we expect significant mod¬ 
ifications of the density shift. 

Sliding Superlattice: We consider the third case where 
the interrogation is done by a pair of counter-propagating 
beams close to resonance with the clock transition and 
with a global phase difference T(t) which can be con¬ 
trolled in time (see Fig. [^a)). The non-interacting 
Hamiltonian, written using a Wannier orbital basis along 
the lattice direction, is 


n,j,a 


Hl — f Jn + l + H.C. 


2 


(^n cos(T(i) - j4)) + H.c. ) . (5) 


As T(t) is changed, the clock laser standing wave “slides” 
with respect to the optical lattice. This allows for the 
minimum realization of a topological pump when T (t) is 
adiabatically varied from 0 —> 27r [231I27H3Q] . In the weak 
tunneling limit the quantized nature of particle transport 
can be directly linked to spatially isolated tunneling res¬ 
onances [23]. We now show how those resonances can 
be spectroscopically measured in OLCs. 
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Let us first consider the case = S = 0 and set 
4> = Ttt/G, relevant for the ^^Sr system. We write 
T = (27r5 +AT)/12, with 0 < AT < 2tt and s an integer. 
Under these conditions, the localized dressed eigenstates 
are spin-polarized along alternating between neigh¬ 
boring sites except for “defects” at = 6r + 3 + 5 

(r an integer) when cos(T — AT/12 — = 0 and the 

ground states at and jd(r) +1 point along the same 
direction. Since tunneling preserves polarization, it is 
suppressed when Jn <C Un due to the energy offset ^ ft 
between neighboring sites. The one exception is the case 
AT = TT where the defect pair jdi^) and jd{'^) + 1 is 
resonantly tunnel-coupled. Quantized transport occurs 
when AT is slowly varied across the resonance. Instead 
of adiabatic transport we propose to spectroscopically re¬ 
solve the resonance using a Ramsey-type protocol. Here, 
atoms initially prepared in ^ at |(5| ^ Uq are adiabat- 
ically transferred to the ground dressed state by slowly 
turning off at T = {27rs + ATp)/12 with AT^ ^ tt. 
Then T is quenched so resonant tunneling is allowed, 
5 ^ s + 1 and AT = tt, and the system evolves for a 
time r. Following this evolution tunneling is turned off 
and the phase switched back, s-\-l ^ s and AT — ATp. 
Those atoms which have tunneled at the resonant sites 
are now in an excited state of the local dressed basis. 
These excitations can be measured by adiabatically con¬ 
verting the dressed excitations to “bare” e excitations by 
adiabatically ramping on (5 ^ Uq; this leads to a mea¬ 
surable excited state population ne{t). In Fig. gc) we 
show these tunneling resonances are clearly observable 
even at finite temperature. Moreover, since resonances 
are spatially well-separated (every six sites), they can be 
resolved with low-resolution imaging. 

Interactions modify the transport dynamics in our pro¬ 
tocol. The sliding superlattice simplifies the treatment of 
interactions by isolating resonant site pairs {jd^jd + !}• 
We consider the case where at most two atoms occupy 
the resonant sites, a condition which can be achieved by 
decreasing the atomic density with a large-volume dipole 
trap m- The two particle states can be classified in 
terms of the atoms’ spin polarization along x in four sec¬ 
tors: three symmetric ones (triplets) with total x polar¬ 
ization Mx = +1,—1,0 and a = 0 singlet. Within 
each sector, the dressed states’ interaction parameters 
are given by [38] U+i = U_i = (U^^ + + 2U®^)/4, 

Vo = {V39 + V^^)/2 and Uq = (see Fig. gb)). 
Since the p-wave parameters are not SU(2) symmetric, 
i.e. bee bgg 7 ^ 6e^, the triplet sectors are coupled. 
However, in the weakly interacting limit U, the 

triplets are separated by energy gaps ^ ft and coupling 
between the sectors can be neglected. The singlet sector 
is always decoupled from the triplets. 

Within each sector interactions modify the dynamics 
making it sensitive to temperature and density. The 
modifications can be extracted by performing measure¬ 
ments of the atom number-normalized excitation frac¬ 
tion for different densities and then differentiating the 
high-density and low-density results [38] . This procedure 


(a) 


o> 








(b) Uo Fo F+i F_i 



FIG. 4: (a) An additional counter-propagating probing beam 
with a differential phase T generates a sliding superlattice 
potential, shown for <f> = Ttt/G, corresponding to an ®^Sr 
OLC. For weak tunneling J ft transport is energetically 
suppressed except at resonant defect points (circled), (b) 
Two-particle interaction sectors classified by the total polar¬ 
ization Mx and spatial symmetry of the dressed states. An 
oval (figure-eight) denotes a symmetric (antisymmetric) spa¬ 
tial wavefunction. (c) Dynamics of a single particle at the 
tunneling resonance for two temperatures (red solid and blue 
dashed) and an off-resonant site (black dotted) for Jo = 8Hz, 
Uo = IkHz, 0 = Ttt/G. (d) Normalized differential excita¬ 
tion extracted as explained in the text with the interaction 
parameters of Ref. m- The Mx = 0 components (solid lines) 
involve the s-wave sector, and so display a strong dependence 
of contrast on temperature, while the M±i sectors (dashed 
lines) experience only weaker single-particle thermal dephas¬ 
ing. 


removes the single-particle contribution and is particu¬ 
larly suitable for characterizing the role of interactions 
in clock experiments [Tain]. The normalized differen¬ 
tial excitation is shown in Fig. id)- For the adiabatic 
dressed state preparation, all four manifolds, and hence 
all interaction parameters 5 contribute to the dy¬ 
namics. A filtering protocol that uses the ee losses can 
be used to separate the dynamics of the various sectors. 
For example, by transferring all atoms to the e state and 
holding before the adiabatic ground state preparation, 
the doubly occupied = ±1 triplet sectors will be re¬ 
moved and only the Mx = 0 singlet and triplet remain 
and contribute (here the ground dressed states have one 
atom at jd and jd + !)• As shown in Fig. [^the dynam¬ 
ics of the U±,o and Uq sectors can be distinguished by 
the different scaling of the p and 5-wave interaction pa¬ 
rameters with temperature Tr, V±p (Tr) ^ const, and 
Uq {Tr) [18] (See FigJ^d)). By comparing these 

dynamics to that without the holding time, information 
about the Mx = ±1 dynamics can be extracted. In gen¬ 
eral 5-wave interactions dominate the normalized differ¬ 
ential contrast, with p-wave contributions relevant only 
at hotter temperatures. 

Summary We have described approaches to implement 
and probe SOC in OLCs. We showed the proposed proto¬ 
cols provide clean signatures of SOC in interacting many- 
body systems. While they work at the temperatures 
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achievable in current OLCs, their sensitivity and applica¬ 
bility are expected to significantly improve when operat¬ 
ing the clock in the quantum degenerate regime. More¬ 
over, if nuclear spins are included, they can open a win¬ 
dow for the investigation of SOC with SU(N)-symmetric 
collisions [43] , 

Acknowledgments We thank Ed Marti and James 
Thompson for comments on the manuscript, and Leonid 


Isaev and the JILA Sr clock team for discussions. This 
work was supported by the NSF (PIF-1211914 and PFC- 
1125844), AFOSR, AFOSR-MURI, NIST and ARO, the 
EPSRC Grant EP/K030094/1, and by the JILA Visit¬ 
ing Fellows Program. MLW thanks the NRC postdoc¬ 
toral fellowship program for support. AK was supported 
by the Department of Defense through the NDSEG pro¬ 
gram. 


[1] V. Galitski and 1. B. Spielman, Spin-orbit coupling in 
quantum gases, Nature 494, 49-54 (2013). 

[2] N. Cooper, Quantum Hall States of Ultracold Gases , In 
Many-Body Physics with Ultracold Gases: Lecture Notes 
of the Les Houches Summer School: Volume 94, 28 June 
- 23 July 2010, Salomon, C., Shlyapnikov, G. V., and 
Cugliandolo, L. F., editors. Oxford University Press, Ox¬ 
ford (2013). 

[3] J. Dalibard, F. Gerbier, G. Juzeliunas, and P. Ohberg, 

Colloquium: Artificial gauge potentials for neutral 

atoms. Rev. Mod. Phys. 83, 1523-1543 (2011). 

[4] N. Goldman, G. Juzeliunas, P. Ohberg, and 1. B. Spiel- 
man, Light-induced gauge fields for ultracold atoms. 
Rep. Prog. Phys. 77 126401 ( 2014). 

[5] H. Zhai, Degenerate quantum gases with spin-orbit cou¬ 
pling: a review Rep. Prog. Phys. 78, 026001 (2015). 

[6] B. K. Stuhl, H.-I Lu, L. M. Aycock, D. Genkina, and 

I. B. Spielman, Visualizing edge states with an atomic 
Bose gas in the quantum Hall regime, arxiv: 1502.02496 
(2015). 

[7] M. Mancini, G. Pagano, G. Cappellini, L. Livi, M. Rider, 

J. Catani, C. Sias, P. Zoller, M. Inguscio, M. Dalmonte, 
and L. Fallani, Observation of chiral edge states with neu¬ 
tral fermions in synthetic Hall ribbons, arXiv: 1502.02495 
(2015). 

[8] X. Cui, B. Lian, T.L Ho, B. Lev, and H. Zhai, Synthetic 
gauge field with highly magnetic lanthanide atoms Phys. 
Rev. A, 88, 011601 (2013). 

[9] Andrew D. Ludlow, Martin M. Boyd, Jun Ye, E. Peik, 
and P. O. Schmidt, Optical atomic clocks. Rev. Mod. 
Phys. 87, 637 (2015). 

[10] M. M. Boyd, T. Zelevinsky, A. D. Ludlow, S. M. Fore¬ 
man, S. Blatt, T. Ido, and J. Ye, Science 314, 1430 
(2006). 

[11] B. J. Bloom, et al., Nature 506, 71 (2014). 

[12] T. L. Nicholson, et al, Nat. Commun. 6, 6896 (2015). 

[13] N. Hinkley, et al, Science 341,1215 (2013). 

[14] U. Ichiro et al., Cryogenic optical lattice clocks, Nat. Pho¬ 
ton., 9, 185-189 (2015). 

[15] M. J. Martin, M. Bishof, M. D. Swallows, X. Zhang, C. 
Benko, J. von-Stecher, A. V. Gorshkov, A. M. Rey, and 
Jun Ye, A quantum many-body spin system in an optical 
lattice clock. Science 341, 632 - 636 (2013). 

[16] N. D. Lemke, J. von Stecher, J. A. Sherman, A. M. Rey, 
C. W. Oates, and A. D. Ludlow, p-Wave cold collisions 
in an optical lattice clock, Phys. Rev. Lett. 107, 103902 
( 2011 ). 

[17] X. Zhang, M. Bishof, S. L. Bromley, C. V. Kraus, M. S. 
Safronova, P. Zoller, A. M. Rey, and J. Ye, Spectroscopic 
observation of SU(N)-symmetric interactions in Sr orbital 
magnetism. Science 345, 1467 (2014). 


[18] A. M. Rey, A.V. Gorshkov, C.V. Kraus, M. J. Martin, 

M. Bishof, M. D. Swallows, X. Zhang, C. Benko, J. Ye, 

N. D. Lemke, and A.D. Ludlow, Probing many-body in¬ 
teractions in an optical lattice clock. Annals Phys. 340, 
311 - 351 (2014). 

[19] M. Bishof et a/.. Inelastic Collisions in Optically Trapped 
Ultracold Metastable Ytterbium, Phys. Rev. A 84, 
052716 (2011). 

[20] A. D. Ludlow, N. D. Lemke, J. A. Sherman, C. W. 
Oates, G. Qumner, J. von Stecher, and A. M. Rey, Cold- 
collision-shift cancellation and inelastic scattering in a 
Yb optical lattice clock, Phys. Rev. A 84, 052724 (2011). 

[21] R. A. Williams et al. Synthetic Partial Waves in Ultracold 
Atomic Collisions, Science335, 314-317 (2012). 

[22] Z. Fu et a/.Production of Feshbach molecules induced by 
spin-orbit coupling in Fermi gases. Nature Physics 10, 
no (2014). 

[23] N. R. Cooper and A. M. Rey, Adiabatic Control of 
Atomic Dressed States for Transport and Sensing, Phys. 
Rev. A 92, 021401(R) (2015). 

[24] Michael Lohse, Christian Schweizer, Oded Zilberberg, 
Monika Aidelsburger, and Immanuel Bloch, A Thouless 
Quantum Pump with Ultracold Bosonic Atoms in an Op¬ 
tical Superlattice, arXiv: 1507.02225 (2015). 

[25] Shuta Nakajima, Takafumi Tomita, Shintaro Taie, Tomo- 
hiro Ichinose, Hideki Ozawa, Lei Wang, Matthias Troyer, 
and Yoshiro Takahashi, Topological Thouless Pumping of 
Ultracold Fermions, arXiv: 1507.02223 (2015). 

[26] Hsin-I Lu et al. Geometrical pumping with a Bose- 
Einstein condensate, arXiv: 1508.04480 (2015). 

[27] D. J. Thouless, Quantization of particle transport, Phys. 
Rev. B 27, 6083 (1983). 

[28] Lei Wang, Matthias Troyer, and Xi Dai, Topological 
Charge Pumping in a One-Dimensional Optical Lattice, 
Phys. Rev. Lett. Ill, 026802 (2013). 

[29] Feng Mei et al, Topological insulator and particle pump¬ 
ing in a one-dimensional shaken optical lattice, Phys. 
Rev. A 90, 063638 (2014). 

[30] Ran Wei and Erich J. Mueller, Anomalous charge pump¬ 
ing in a one-dimensional optical superlattice, Phys. Rev. 
A 92, 013609 (2015). 

[31] F. Scazza, et al. Observation of two-orbital spin-exchange 
interactions with ultracold SU(N)-symmetric fermions 
Nature Phys. 10, 779-784 (2014). 

[32] C. Zhang, S. Tewari, R. M. Lutchyn, and S. Das Sarma, 
PxPipy Superfluid from s-Wave Interactions of Fermionic 
Cold Atoms, Phys. Rev. Lett. 101 , 160401 (2008). ? 

[33] B. Julia-Dfaz, T. GraB, O. Dutta, D.E. Chang, and M. 
Lewenstein, Engineering p-wave interactions in ultracold 
atoms using nanoplasmonic traps, Nat. Comm. 4, 2046 
(2013). 



6 


[34] F. Cappelinni, et al. Direct observation of coherent inter¬ 
orbital spin-exchange dynamics Phys. Rev. Lett. 113, 
120402 (2014) 

[35] Michael L. Wall, Kaden R. A. Hazzard, and Ana Maria 
Rey, Effective many-body parameters for atoms in non- 
separable Gaussian optical potentials, Phys. Rev. A 92, 
013610 (2015). 

[36] A. Cell, P. Massignan, J. Ruseckas, N. Goldman, LB. 
Spielman, G. Juzeliunas, and M. Lewenstein, Synthetic 
gauge fields in synthetic dimensions, Phys. Rev. Lett. 
112, 043001 (2014). 

[37] Dario Hiigel and Belen Paredes, Ghiral Ladders and the 
Edges of Ghern Insulators, Phys. Rev. A 89, 023619 
(2014). 

[38] See the Supplementary Material at [xxx] for a derivation 
of the density shift in Ramsey spectroscopy, the extrac¬ 
tion of the chiral Bloch vector angle in Rabi spectroscopy, 
and the interaction-induced contrast decay with the slid¬ 
ing clock superlattice. 

[39] Y. Hatsugai, Topological aspects of the quantum Hall 
effect, J. Phys.: Gondens. Matter 9, 2507 (1997). 

[40] L. Tarruell, et al.Greating, moving and merging Dirac 
points with a Eermi gas in a tunable honeycomb lattice. 
Nature 483, 302 (2012). 

[41] M. At ala, et ah. Direct measurement of the Zak phase in 
topological Bloch bands, Nat. Phys. 9, 795 (2013) 

[42] L. Duca, et ah. An Aharonov-Bohm interferometer for 
determining Bloch band topology. Science 347, 288 
(2015) 

[43] M. Gazalilla and A. M. Rey, Ultracold Eermi Gases 
with Emergent SU(N) Symmetry, Rep. Prog. Phys. 77, 
124401 (2014). 



7 


Supplemental Material for “Exploring synthetic spin-orbit coupling in a 

thermal optical lattice clock” 

Interaction Hamiltonian with the spin model assumption The Hamiltonian governing fermionic AEAs in an optical 
lattice clock may be written: 






A 


(r) — h5 JI 


dripl (r) xl)e (r) , 


m 

TT 

Hl = -^ 


j drV)| (r) tAe (r) (r) (r) + y] j drW (r), (w t/)^(r)y^(r) y 

al3 

J dr [t/.t (r) (r) 


+ H.C. 


(51) 

(52) 

{S3) 

(S4) 


where m is the atomic mass, (i*) is a fermionic field operator for state a G e}, W A (r), H (r) = (VA {t))B (r) — 


A (r) (VH (r)) is the Wronskian, and (5 = cj/ — c^o is the difference between the clock laser frequency uji and atomic 
frequency ujq in the rotating frame of the laser. The Hamiltonian Hq contains the kinetic energy and trapping 
potential fAxt (i*)? Hj contains the effects of s-wave interactions with scattering length a~g and p-wave interactions 

with interaction volumes 6^^ between nuclear spin-polarized AEAs, and Hl describes the coupling of internal levels to 
the clock laser. We assume that the optical lattice is deep enough that we can neglect interactions occurring between 
different lattice sites, and so consider only interactions between particles occupying the same lattice site. Eurther, 
we enact the spin model approximation [18], which neglects all interactions which do not preserve the single-particle 
transverse mode occupations during a two-body collision. The spin model Hamiltonian hence keeps only direct terms, 
in which the two colliding particles remain in their same transverse motional quantum states, and exchange processes 
where the transverse quantum numbers of the two colliding particles are exchanged. The spin model Hamiltonian is 
valid when interactions are smaller than the spacing between transverse modes, and also when anharmonicity of the 
potential prevents collisional exchange of mode energy between transverse dimensions. 

To derive a Hubbard model amenable for direct calculation, we expand the field operators in a basis of single-particle 
eigenfunctions. In order to facilitate thermal averages, we take this set of functions to be the eigenfunctions of 

(r), which are indexed in terms of a quasimomentum q and a transverse mode index n. Enacting this expansion 
with the approximations of the last paragraph, we find 


= T. E [( 2 -''"."=)«.,..)+<.=>)< 


g+Ag,ni^^,g'-Ag,n2^/5»g'>n2^«,g,ni 


a/3 {ni,n 2 } qq' 


cx(3 


+ 2(1 (^nin2) (y 


-u\ 


cxf3 


{ni,n2} ^{ni,n2} 


}) 




-Ag,ni^a,g+Ag,n2^/559',n2<^a,g,ni 


(S5) 


Here, L is the number of lattice sites, the sum over {ni, n 2 } means the sum over distinct, ordered pairs of transverse 
modes, g, g', and Ag are quasimomenta in the first Brillouin zone (BZ), greek letters denote internal electronic states 
G {g, e}, and the interaction matrix elements may be written as 


U' 


{ni,n2} ~ dsS, 




m 


^nin2n2ni , ’^{ni,n2} 




al3 _ 


127rfi^ 


m 


Ir P 

^a/3^nin2n2ni 


(S6) 


by defining the integrals 


Sn' 


= J (r) V’/.n' (r) V’i,n2 (r) z/’i.m (r) , = J 


= / drW 


(r), V’/.n' (r)] W (r), (r)] . 

(S7) 


In the latter two expressions, the wavefunctions are localized Wannier-type functions obtained from 
the quasimomentum-indexed eigensolutions (Bloch functions) of the full 3D potential as (r) = 




















these functions as 


We can also define the mode-dependent tunneling and Rabi frequency in terms of 




J„ = - y drtpj^a (r) 

fln = ^ J dri^j,n (r) (r) 


V’i+i.n (r) , 


(58) 

(59) 


First-order density shift in Ramsey speetroseopy In this section we compute the density shift for Ramsey spec¬ 
troscopy to first order in interactions using an interaction picture perturbation series. In particular, we write the propa¬ 
gator during the dark time evolution in terms of a truncated Dyson series i . 

At the end of the dark time, a spin-rotation pulse of area O 2 is applied, and then the total is measured. These 
two operations can be combined as the measurement of the operator Sq = Yhnq 
first-order result of the Ramsey sequence is 


cos OSq^ — sin 


Hence, the 


(5^ (r)) = (V;(0i) (ei))] , (SIO) 

L Jo 

where \ fj (6>i)) is the state resulting from applying the first Ramsey pulse to an initial state \i). We will parameterize 
the initial state in terms of products of single-particle eigenstates \i) = Y[f=i where |vac.) is the vacuum, 

N the number of particles, and the initial distinct set of populated quasimomenta and transverse modes, so 

that 


/ 0 0 \ 

l^(^'i)) = n (^cosyaj +isinya^^„,j . 


(Sll) 


From this, we find that the non-interacting dynamics are given as 

N 




i=l ■- 


COS 0i cos O 2 sin 0i sin O 2 


cos {6t - AEm {qi,4>) t) 


where (q) is the single-particle dispersion and 

AEn {q, </>) = En{q+ <!))- E^ (q) 


(S12) 


(S13) 


is the difference in single-particle dispersion of the e and g states. 

Writing the non-interacting result as {S^ (r)) = A (r) cos (Sr) + B (r) sin (Sr) + const, we can extract the density 

shift as An = 2 ;^ai‘ctan and the normalized contrast decay as C (r) = . For a single particle 

with momentum q and transverse mode n, the density shift is AEn {q^ /27r. Assuming short times compared to the 

tunneling bandwidth, the density shifts for all particles add as An = AEm {qi^ 0) I^ttN + O (t^). The contrast 

decay at short times is given as C (r) = 1 - {qi, </>)) + w (fti </>) -|- O (t^). This 

single-particle contribution will dominate over the ^ U‘^r^ decay of the contrast due to interactions at short times 
and in the regime of weak interactions compared to the tunneling. 

The contribution to the dynamics at the lowest order in interactions is given as 


Jo 

^ I sin $1 sin O 2 sin 


(S14) 


{Pl,P2} 


J^Enp {qpi , 4>) + AEnp {qp2 ? 


(C'{n,,n,J-X{n,,n,jCOS0l)cOs(^ 
COS 0i 

- On.. ^ 1 COS 


(2 np2 ) C|np^ np2 } 


J^Enp {qp^ , (f) AEnp {qp2 5 0) 


/ AEnp {qp^^cj)) AEnp {Qp2^^) 


sine I 


/ AEnp J^Enp (gp 2 5 0 ) 


]]} 
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where the mode-dependent spin model parameters are 


C. 


{nin2} 


X{nin2} 


C{nin2} 


T/ee _ ygg 
^{nin2} ^{nin2} 


vr X + . - 2vr X 

{nin2} {nin2}{nin2} 


yeg _ jjge 
{nin2} |nin2j 


(515) 

(516) 

(517) 


Taking a series in the tunneling bandwidth, the zeroth order term for a given pair of particles is 
(C — xcos^i) sin^i sin ^2 sin(5r, which has been obtained in previous works [15]. The thermally averaged density 
shift in the high-temperature limit will be given in the next section. 

High-temperature thermal average of density shift Here, we derive the high-temperature (compared to the tunneling 
bandwidth) thermal average of the first-order many-body dynamics derived in the last section. First, we expand the 
above result to third order in time, finding 


1 ^ 

(S' (T)) = 5 E (cos (St) sin sin O 2 — cos 0i cos O 2 ) 


(SI 8 ) 


i=l 

+ E 

{Pl,P2} 

- E 

{Pl,P2} 
N 


sin 0i sin O 2 


sin 0i sin O 2 
L ~2 


1 ^ 

^ COS^i) silK^T + 2 E 

i=l 

,np2 } ~ ,np2 } ^ ’ ^)) 


- j ^ (qi, 0)^ COS (Sr) sin 6 >i sin 6 > 2 r^ 




- E 

{Pl,P2} 


sin 0i sin O 2 r 


4 ^ (^{npi,np2} ^{npi,np2} + ^^^^2 {Qp2^^) ) 


L 24 L 

— (2 — Srip^ np2 ) ,np2 } ^^P^ ’ ~ ^^^ P 2 ^^P‘^ ’ 

1 ^ 

- (gi,(/))^sin((5r)sin6>i sin6>2T^| . 




To perform thermal averages, we now enact the tight-binding approximation, in which E^i (q) = E^i — 2 Jn cos (g), and 
compute the partition function as 


Z = '^e ~ E^ 

nn VI J —1 


cos(7rg) 


= 2 ^ Jo (2/?J„) = 2 ^ e-^®" [l + O {/3^J^)] , 


(519) 

(5 20 ) 


where (x) is the modified Bessel function of order n and [3 the inverse radial temperature. Hence, for temperatures 
sufficiently high compared to the bandwidth Jn, the partition function is just a constant times the partition function 
of the transverse modes. Now, 

(AE„ (</, 0))t« = I E J„ (1 - cos 0) h (2/?J„) 

n 

= TEe“^^"2/3J„(l-cos(/;)) + C)(/?2j2) , 

„ 

where Z„ = is the partition function of the transverse modes. This result gives that terms which are linear 

in AJln {q, <p) for a specific q, including terms like AEm {qi,4>) ( 92 , (f>), vanish at least as fast as O (/?) at high 


(5 21 ) 

(5 22 ) 
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temperatures. In contrast, we find 


{AEl (q, 4 >))t^ = (2/3 J„) + I 2 (2/3 Jn) (1 - cos <P) 

= (1 - cos c/)) + 0 {P^Jl) , 


. 2 

Sin - 


(S23) 


and so the thermal average of {q, (j)) is a constant to lowest order in a high-temperature expansion. 
Using these results, we find that the leading order approximation in a high-temperature series expansion is 

{{^e2))TR = (sm< 9 i sin<92 cos (( 5 t) - cos( 9 i cos6>2) H- — -^sin6>i sin6>2 sin (( 5 r) {{C)tr - {x)tr cos6>i) 


2 

2N{J‘^)tr sin 6 >i sin 6>2 cos {5t) sin^ ^ 


2L 

2^2 (JV„r 3 iV (TV-!).(/. 


sin - [6 ((C')rH - (x)rK cos6»i) - 2{C)t^ cos6»i] 

(S24) 


where we have neglected interactions between pairs of particles with the same transverse mode for simplicity, {•)tr de¬ 
notes a thermal average with respect to the transverse modes, and we have assumed that, e.g. {J‘^C)tr ~ {J‘^)tr {C)tr^ 
which is valid when the tunneling is only weakly temperature dependent. From this, we find the density shift 


Az/ = Az/( 


1 I 4 , 2 ^ _2 ■ 2 '/> (C)r„cos0i 

1+3^ )t,t ,m ■ 


(S25) 


where Az/q is the density shift in the absence of tunneling. 

Lineshape and extraction of O^q from momentum-resolved Rabi spectroscopy 

As described in the main text, SOC introduces substructure in the Rabi lineshape within a window of width 
8 Jn sin ^ of the carrier frequency. At finite temperature, many transverse modes are populated and hence the 

dependence of Un and Jn on n could broaden the line and destroy this substructure. However, a direct simulation of 
the Rabi lineshape using the potential Eq. 0 demonstrates that the features of the ideal, zero-temperature lineshape 
are captured up to a temperature of 3/jK as shown in Fig. [^a). Thanks to the clock’s sub-Hz resolution it will be 
possible to resolve the lineshape for a wide range of lattice parameters. 

We now turn to the extraction of the chiral Bloch vector angle O^q with Rabi spectroscopy. The eigenstates of 
Eq. 0 may be written as 


with energies 


IV’n.g,-) = cos^lgqn) +sm^\eqn) 

\'ipn,q,+) = - sin ^ \gqn) + cos ^ \eqn), 


p _Eniq) + En(q + (p) + S 


(526) 

(527) 

(528) 


and magnetization 


/ / I I / \ _ ^ 

('0n,g,± |5^n,gl'^n,g,±) — ^ 


(S29) 


Using the above, the expected Rabi dynamics in the ideal case for an arbitrary initial state a\grvq) -\- /3|eng) is found 
to be 


(a,/3,t) = i ^|/3|^ - [cos^ 6ng + cos ( 2 B„,ji) sin^ + TZ (a*/3) cos 6ng sin 6nq (cos( 2 B„gT) - 1) 

+1 (a* 13) sin Onq sin ( 2 B„qt) . (S30) 

The protocol extracting the chiral Bloch vector angle O^yq involves three sequences (Fig. I^b)). All sequences start 
by preparing the atoms in g. In sequence I, a narrow 7 r-pulse about the x axis is applied at the detuning 5'^ associated 
with the q^ resonance. The Rabi frequency of the pulse, Oq, should be weak enough to guarantee that only atoms 
within a narrow window centered around q^ are transferred to e. Atoms with qi^ q^ are off-resonant and remain in g. 
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FIG. 5: (Color online) (a) Thermally averaged Rabi lineshape fully accounting for trap non-separability. SOC-induced peaks 
are visible even at T = 3/iK. (b) Three pulse sequences used to extract Ooq in Fig. (2) of the main text. 


Next the detuning and Rabi frequency are quenched to the desired values Qq and and Rabi oscillations are recorded. 
Sequence II is identical, except that the initial pulse is about —x. Sequence III uses no initial pulse, but is otherwise 
the same as sequence I. The dynamics in I contains information about Oq *, but this information will be buried in the 
signal of the q ^ Qq atoms. Subtracting the dynamics of I and III isolates the dynamics of Qq atoms. Sequence II 
is required because due to the mode dependence of fin, a 7r-pulse is not experienced by all atoms populated at the 
transverse modes at high temperature. The dynamics of the atoms with Qq that remain in g is cancelled by taking 
the average between sequences I and II, which effectively deals with only the atoms transferred to e. The dynamics 

for the sequences I, II, and III given in the main text are i/n, , and (1, 0, t), 

respectively, where /n ~ 1 is the excitation fraction resulting from the initialization 7r-pulse. This gives the subtracted 
signal 


S'" - /n> -ifn, t) + 5"" - /n> '^fn, *) 


- (1, 0, t) = (cos^ Onq + COS {2\Bng\t) siii^ 6>„,) 


(S31) 


For a single frequency |5ng|, and O^q can be inferred from these oscillations as 2tan^^ng = (max — min)/mean. 
In the multiple-frequency case this is no long er a strict equality, but Fig. I^b) shows that determining O^q in this 


fashion from the thermal average of Eq. (S31) is nevertheless robust. In addition, the validity of this approach can 
be quantitatively estimated by the equality (max + min)/mean = 2, which employs the same assumptions. 

Interacting dressed states with a sliding clock superlattice Here, we consider the dynamics of two atoms in neighboring 
sites jd(r) and + 1 which are tunnel-coupled for a sliding lattice clock phase AT = tt. We will take these two 
atoms to have different transverse mode indices n and m. Interactions in the s-wave channel occur between these 
atoms when they are in an antisymmetric electronic state |A^) = {\ge) — \eg)) /\/2 with a symmetric spatial wave 
function |^nm) = with the subscript j denoting the lattice site index. Instead, p-waye interactions 

occur when the atoms are in a symmetric electronic state |A+) = {\ge) + \eg)) /\/2 and an antisymmetric spatial 
wave function |^nm) = purposes of describing the system dynamics at AT = tt, where 

I ~^) = (1^) + k)) /v^ is the single-well electronic ground state and | = {\g) — |e)) /\/2 the excited state, it is 

useful to employ the states 


t, 1 : | 2 ^^, 0 ), | 0 , 2 ^^), I ^n, ^m), I ^m, ^n) , 
t,— 1:|2 , 0 ),| 0,2 ), I ^n, ^m), I ^m, ^n) 


(S32) 


i,0: |2±,0),|0,2±), 
s,0: |2T,0),|0,2^), 


^115 m) I n) | ni ^m) | ^ mi ^n) 

■ni ~^m) + I ^mi ~^n) 


72 

I ■^m) + I ■^n) 


72 


72 


|2++) = |^^)|^-m), 

|2±) ^ + I ^ 

|2^) = f 


72 

-)-| 


72 


I ^ nm/ 




(533) 

(534) 

(535) 

(536) 


where 
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In each four-state sector the Hamiltonian may be written as 








0 

J 

J 


0 

-J 

-J 


J J 
-J -J 


(S37) 


where ja = Us,o = Uq, I is the identity operator, and we have set ^ and Jn = Jm = J 

for simplicity. The s, 0 singlet sector is rigorously decoupled from the others, but there is mixing between the triplet 
sectors proportional to the spin model parameters y and C. These couplings are neglected due to the large single¬ 
particle energy difference ^ between coupled manifolds with the assumed separation of energy scales ^ J, V, [/. 
The Hamiltonians Eq. (S37) may be readily diagonalized, leading to the eigenvalues 2Ma,0,2Ma,ll + t/, 2Ma,0 


^ At high temperatures compared to the tunneling, and for the adiabatic preparation procedure 

described in the main text, the observed dynamics will be an equal weight superposition of the dynamics from the 
initial states |2++,0), |0,2 ), | and | The dynamics of the number of excitations following 

adiabatic conversion back to the “bare” gje basis for these states are 


( I W\ n sin ^ sin (iy4J2~+V^74) 

|2++,0),|0,2 ) ^ne(t) =l-cos —cos fy4J2 + —j- ^IQJ2 + ]/2 - 


Vot 


—^11) ■^m), I —^m) ■^n) “^ IT-e (t) —1 — - [ COS COS I t\ 4J2 H ^ ] + COS —— COS I t\ 4J2 H ^ 




Uot 




(S39) 


Vo sin ^ sin ViWVFV) Uq sin ^ sin (t y4J2 + f/2/4) ' 


yi6 j2 + q2 


yi6j2 +1/2 


For small interactions compared to tunneling, we can expand the result to lowest order in interactions, and find the 
thermally averaged dynamics 


O • 2 , {Ui + V? + 2FV t (sin {2Jt) - 2Jt cos {2Jt)) 

Ue [t) = 2 Sin [Jt) - ^^- — -. (S40) 

In addition to the above two-particle dynamics, there will be contributions from experimental realizations in which 
the resonant double well has only a single particle. Here, the dynamics is rig (t) = sin^ (Jt). Writing pi as the 
probability of a single particle in a given double well and p 2 as the probability of having two particles, the small- 
interaction dynamics averaged over all realizations is hence 


ne (t) = (pi + 2p2) Sin^ {Jt) - p2 






2E]^) t (sin (2 Jt) — 2 Jt cos (2 Jt)) 


64J 


(S4I) 


Noting that pi 2p2 = N the average number, we can subtract the measurements of Ug (t) /N for two different 
densities with single- and double-occupancy probabilities (pi,P2) and {p'i^p' 2 ) and numbers N and N'^ respectively, 
to find 


A[ne {t)/N] = -{?^ 


P 2 ^ {^0 + V + t (sin {2Jt) - 2Jt cos {2Jt)) 
N') 64J 


(S42) 


In this way we separate the single-particle contrast decay due to a thermal spread in J„ from the contrast decay due 
to interactions. 




















